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ABSTRACT 

Analytical,  numerical  and  experimental  results  are  compared  in  a 
systematic  investigation  of  the  effects  of  projectile  shape  on  thin  plate 
perforation.  Ihe  analysis  of  five  low  velocity  impacts  involving  optimal-, 
ogival-,  Russian-,  conical-  and  cylindrical-head  projectiles  comprise  the 
study.  A  three-way  comparison  of  simple  perforation  theory  predictions, 
HELP  numerical  solutions  and  available  experimental  data  is  made  for  the 
five  problems. 

Ihe  HELP  solutions  are  in  excellent  agreement  with  the  experimental 
data.  Although  the  simple  perforation  formulas  are  sensitive  to  pro- 

i 

jectile  shape,  they  do  not  provide  adequate  residual  energy  predictions. 
Further  numerical  studies  of  this  type  should  provide  additional  insight 
into  the  penetration  process  and  eventually  lead  to  improved  perforation 
formulas . 
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I.  INTRODUCTION 

A  few  basic  considerations  suffice  to  indicate  the  moti¬ 
vation  for  the  present  numerical  analysis  of  projectile  shape 
effects.  The  efficiency  with  which  a  kinetic  energy  penetrator 
is  able  to  defeat  a  particular  target  is  dependent  on  the  con¬ 
figurational  and  material  parameters  as  well  as  the  mode  of  plate 
failure.  Existing  perforation  formulas  cannot  accurately  pre¬ 
dict  this  dependency.  Furthermore,  the  application  of  these 
formulas  requires  that  the  plate  failure  mode  be  assumed  in 
advance.  Numerical  techniques,  however,  can  be  used  as  a  tool 
to  investigate  the  complicated  relationships  between  material 
and  configurational  parameters,  material  flow  history,  failure 
criteria  and  plate  failure  mode.  Such  techniques  should  there¬ 
fore  be  useful  in  predicting  penetration  efficiency  and  provide 
useful  information  to  designers  of  both  weapons  and  armor  as 
well  as  to  vulnerability  analysts. 

Penetration  efficiency  is  a  difficult  parameter  to  de¬ 
fine  in  general;  however,  for  the  case  of  thin  plate  perforation, 
it  can  be  defined  as  the  ratio  of  residual  projectile  kinetic 
energy  to  initial  energy.  Even  with  such  a  precise  definition, 
penetration  efficiency  is  difficult  to  calculate  analytically 
since  it  requires  a  perforation  formula  which  accurately  pre¬ 
dicts  the  residual  projectile  velocity  and  which  is  extremely 
sensitive  to  the  configurational  and  material  parameters  which 
define  a  particular  impact.  Furthermore,  since  failure  modes 
also  depend  on  impact  parameters,  a  particular  thin  plate  per¬ 
foration  formula  can  at  best  be  expected  to  provide  good  pre¬ 
dictions  only  over  a  narrow  range  of  parameters.  For  example, 
an  increase  in  projectile  hardness  should  increase  penetration 
efficiency;  however,  if  the  impact  occurs  in  a  regime  in  which 
projectile  break-up  or  shattering  occurs,  penetration  efficiency 
decreases.  Another  example,  and  one  that  will  be  discussed  in 
the  following  sections  of  this  report  in  more  detail,  involves 
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the  effect  of  projectile  shape.  Thin  plate  failure  modes  can 
vary  from  petalling  to  plugging  as  the  projectile  nose  is 
varied  from  an  extremely  pointed  to  an  extremely  blunt  shape. 
Combined  petalling  and  plugging  modes  can  occur  for  inter¬ 
mediate  shapes.  Since  projectile  shape  is  an  important  para¬ 
meter  to  which  thin  plate  penetration  efficiency  is  more  sen¬ 
sitive  than  can  be  demonstrated  by  existing  perforation  formulas, 
and  since  it  represents  a  parameter  which  when  varied  encompasses 
more  than  one  plate  failure  mode,  it  has  been  chosen  as  the 
subject  of  the  numerical  impact  investigation  reported  here. 

In  Section  II  various  existing  thin  plate  perforation 
theories  are  presented.  In  Section  III  some  recent  work  by 
T.  E.  Fields ^  concerning  the  comparison  of  some  of  the  per¬ 
foration  formula  predictions  with  experimental  data  for  various 
projectile  shapes  is  discussed.  The  inadequacy  of  certain  of 
these  formulas  to  predict  experimentally  observed  residual  veloc¬ 
ities  becomes  quite  apparent.  Numerical  solutions  to  five 
problems  chosen  to  reproduce  selected  impacts  from  Fields'  ex¬ 
perimental  matrix^  are  presented  in  Section  IV.  The  pro¬ 
jectile  shapes,  which  will  be  defined  later,  are  referred  to  as 
optimal-,  ogival-,  Russian-,  cylindrical-  and  concial-head 
projectiles.  The  HELP  computer  code  was  employed  in  these 
solutions.  The  numerical  techniques  are  briefly  described  and 
the  material  model  employed,  including  the  equation  of  state, 
deviatoric  constitutive  relation,  yield  condition  and  failure 
criterion  is  discussed.  The  HELP  solution  to  the  Russian  pro¬ 
jectile  impact  is  presented  in  detail  while  only  selected  numerical 
data  from  the  solutions  of  the  other  four  problems  are  included. 

In  Section  V  the  numerical  results  are  discussed  and  compared 
with  both  the  predictions  from  the  elementary  perforation 
formulas  and  the  experimental  observations. 

Finally,  Sections  VI  and  VII  respectively  present  the 
conclusions  of  the  current  study  (performed  under  the  technical 
direction  of  Dr.  R.  T.  Sedgwick)  and  outline  suggestions  for 
additional  related  research. 
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This  report  is  the  third  under  the  current  contract. 
The  first ^  applies  the  HELP  code  to  two  specific  test  prob 
lems.  The  problems  were  specified  to  provide  a  basis  of  com 
parison  of  the  HELP  approach  with  Lagrangian  methods.  The 
second  report^  ;  is  a  complete  documentation  and  listing  of 
the  HELP  code. 

We  believe  that  the  results  reported  here  are  very 
encouraging  and  point  to  the  desirability  of  additional 
applications  of  the  HELP  code  to  armor  penetration  studies. 
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II.  •  SIMPLE  PERFORATION  FORMULAS 

A  penetration  formula  can  be  defined  as  an  expression, 
theoretically  or  empirically  derived,  relating  the  configura¬ 
tional  and  material  properties  which  characterize  a  particular 
impact  situation  to  those  of  the  residual  state.  For  thick 
target  impacts,  where  the  projectile  is  defeated,  it  is 
customary  to  describe  the  residual  state  by  a  single  para¬ 
meter,  penetration  depth  (the  state  of  the  defeated  pro¬ 
jectile  is  usually  not  of  interest).  For  thin  plate  per¬ 
foration,  the  residual  state  should  include  the  kinetic 
energy,  velocity,  mass  and  shape  of  all  projectile  and  target 
fragments  as  well  as  the  final  target  hole  size.  In  the 
interest  of  simplicity,  however,  most  thin  plate  theories 
are  restricted  to  describing  the  residual  state  by  the 
kinetic  energy  of  the  major  residual  particles,  only.  This 
simplification  is  usually  arrived  at  by  assuming  rigid 
projectile  penetration  in  conjunction  with  an  a  priori 
choice  of  hole  size  and  plate  failure  mode. 

A  brief  summary  of  low  impact  velocity  penetration 
mechanics  is  given  in  Ref.  4.  In  the  following  paragraphs, 
existing  thin  plate  perforation  formulas  pertinent  to  im¬ 
pact  situations  of  interest  in  the  current  study  are  pre¬ 
sented.  Since  the  derivations  of  these  formulas  are  well 
documented  elsewhere^  ^  only  brief  discussions  of  their 
nature  and  applicability  as  well  as  the  physical  assump¬ 
tions  involved  in  their  development  are  presented  with  each 
perforation  formula. 
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2.1  THOMSON  MODEL 

The  perforation  formula  proposed  by  Thomson is 
based  on  the  assumption  that  the  reduction  in  projectile 
kinetic  energy  is  a  result  of  the  work  done  by  the  projec¬ 
tile  in  both  plastic  deformation  and  radial  motion  of 
the  displaced  target  material.  The  projectile  is  assumed  to 
be  rigid  and  the  target  material  is  assumed  to  behave  in  a 
rigid-perfectly  plastic  manner.  In  addition,  the  radial  and 
axial  stresses  are  assumed  small,  relative  to  circumferen¬ 
tial  stresses.  Finally,  it  is  assumed  that  the  hole  size 
equals  the  maximum  projectile  radius,  that  plastic  deforma¬ 
tion  occurs  without  a  change  in  volume,  and  that  energy  loss 
due  to  frictional  forces  at  the  projectile-target  interface 
are  negligible. 

For  the  variables  defined  in  Fig.  1,  the  reduction 
in  projectile  kinetic  energy,  according  to  Thomson's  theory, 
can  be  expressed  as 


W 


irR2h  av 

— y-* 


.h  v2 


^dy 

dx2 


+  2 


/ 


dy 


(1) 


where  a  and  p  are  the  target  yield  strength  and  mass 
density  respectively  and  V  is  the  velocity  of  penetration 
which,  for  simplicity,  is  assumed  to  be  constant  and  equated 
to  the  impact  velocity,  V  .  The  expression  involving  residual 
velocity,  Vr,  is  obtainable  from  the  energy  balance 

„  .  1  Mp  -  \  Mp  V’  C2) 
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Fig,  1- -Projectile  target  configuration  for  describing  the 
geometric  variables  and  parameters  employed  by  the  Thomson 
formula.  The  coordinate  system  is  attached  to  the 

projectile. 
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and  becomes 


irR2h  o 
o 


where  Mp  is  the  projectile  mass. 


(3) 


2.2  ZAID  AND  PAUL  MODEL 

Zaid  and  Paul ^  in  their  model  make  use  of  conserva¬ 
tion  of  axial  momentum  during  rigid  projectile  perforation 
of  a  thin  plate  assumed  to  fail  by  petalling.  Figure  2 
illustrates  the  geometrical  parameters  and  variables  employed. 
Conservation  of  axial  momentum  is  given  by 

MpVo  -  MpV  +  Mt(x')  (4) 

where  V  is  the  instantaneous  projectile  velocity  and  Mt(x) 
is  the  axial  component  of  plate  momentum.  For  pointed  projec¬ 
tiles  the  total  axial  momentum  imparted  to  the  plate,  Mt,  can 
be  expressed  as 


The  equation  for  projectile  velocity  loss  then  becomes 


(6) 
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Fig.  2- -Projectile-target  configuration  for  describing 
the  geometric  variables  and  parameters  employed  by 
the  Zaid  and  Paul  formula.  The  coordinate  system  is 
attached  to  the  projectile. 
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where  the  instantaneous  velocity,  V,  of  Eq.  (5)  has  been 

approximated  by  the  initial  impact  velocity,  V  . 

'0 

2.3  NISHIWAKI  MODEL 

The  perforation  formula  developed  by  Nishiwaki^  is 
arrived  at  by  writing  the  equation  of  motion  of  the  projectile 
acted  upon  by  plate  forces  due  to  a  static  and  dynamic  pressure 
as  well  as  frictional  forces.  For  the  conical  head  projectile 
shown  in  Fig.  3,  the  projectile  equation  of  motion,  assuming 
no  frictional  forces,  can  be  written  as 


M 


dV 
p  <Tt 


♦  pV2sin20 jsinedA' 


(7) 


where  A'  is  the  projectile-plate  contact  area  and  is  an 

empirically  determined  static  pressure  component.  The  ex¬ 
pression  for  the  residual  velocity  for  the  conical  impact 
case  of  Fig.  3  can  be  obtained  by  integrating  Eq.  (7)  and  is 
given  by  Fields ^  as 
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psin2e 
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psin2e 
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2.4  RECHT  AND  IPSON  MODEL 

Recht  and  Ipson^  assume  a  plugging  failure  mode  and 
use  conservation  of  momentum  and  energy  to  express  the  reduc¬ 
tion  in  projectile  velocity  due  to  both  the  acceleration  of 
the  plug  mass  and  shear  resistance  at  the  plug  periphery,  in 
terms  of  the  geometrical  impact  parameters  of  Fig.  4,  the 
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Fig.  3- -Projectile-target  configuration  for  describing 
the  geometric  variables  and  parameters  employed  by 
the  Nishiwaki  formula  for  the  case  of  a  conical 

projectile. 
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Fig.  4*-Projectile-target  configuration  for  describing  the 
geoaetric  paraaeters  eaployed  by  the  Recht  and  Ipson 
foraula  for  the  case  of  a  cylindrical  projectile. 
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expression  for  residual  velocity  is 


where  fl  is  the  ratio  of  target  to  projectile  de*»*ity  and 
V  is  the  minimum  perforation  velocity.  Inherent  in  the 

*  i  1 

development  of  Lq.  (9)  is  the  assumption  that  the  average 
dynamic  shear  stress  acting  on  the  plug  periphery  reaains 
constant.  It  should  also  be  Mentioned  that  the  value  of 
Vxn  must  be  determined  experimentally  and  an  assumption 
regarding  the  plug  diameter  must  be  made. 


2 . 5  HEVDA  MO OIL 

Hcyda1^  has  recently  developed  a  thin  plate  plugging 
model  which  reduces  to  the  theory  of  Kecht  and  Ipson  for  the 
case  of  like  material  impact.  Heyda's  expression  for  residual 
velocity  can  be  written  as 


(10) 


where  the  subscripts  t  and  p  refer  to  the  target  and 
projectile  material  respectively  and  c  is  the  hydrodynamic 
sound  speed.  Heyda  also  gives  an  expression  for  the  minimum 
perforation  velocity,  Vxn,  as 


xn 
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where 


n 


2p,c 
P  P 


ptct * 


°pcp 


> 


K 


n 


and  1  can  be  taken  as  the  length  of  the  projectile  at  the 
tiae  the  initial  pressure  wave  is  attenuated  by  the  rear  free 
surface  rarefaction. 


13 


3SR-348 


III.  COMPARISON  OF  FORMULAS  WITH  EXPERIMENT 

The  following  comparison  of  perforation  formula  predic¬ 
tions  with  experimental  results  is  based  on  the  work  of  T.  E. 
Fields. ^  Fields  investigated  experimentally  the  effect  of 
projectile  shape  on  thin  target  perforation.  Five  different 
shapes  were  analyzed  in  this  study  and  are  referred  to  here 
as  ogival-,  optimal-,  Russian-,  conical-,  and  cylindrical- 
head  projectiles.  The  relative  shapes  are  shown  in  Fig.  5. 

The  nominal  impact  velocity  for  the  set  of  experiments  was 
approximately  S  x  10*  cm/sec.  The  projectile  masses  were  all 
about  20  g  and  both  steel  and  aluminum  targets  were  impacted. 

Three  thicknesses  of  aluminum  plates,  namely  0.159, 
0.317,  and  0.476  cm,  were  impacted  while  only  one  thickness, 
0.317  cm,  was  chosen  for  the  steel  plate  impacts.  Fields 
compared  his  experimental  results  with  the  predictions  of  the 
perforation  formulas  of  Thomson  and  Nishiwaki.  A  sample  com¬ 
parison  for  each  projectile  shape  is  given  in  Table  I  where  it 
is  clearly  evident  that  agreement  is  less  than  acceptable. 

The  bracketed  values  in  Table  I  were  calculated  by  A.  J.  Good. 
The  remainder  of  the  table  is  taken  from  Ref.  1.  it  was  the 
inability  of  the  perforation  formulas  to  predict  the  observed 
residual  velocities  that  lead  to  the  numerical  projectile  shape 
effects  study  discussed  in  Section  IV.  The  lack  of  agreement 
is  not  surprising  in  view  of  the  many  over-simplifying  assump¬ 
tions  required  in  the  development  of  a  simple  perforation 
formula. 
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IV.  NUMERICAL  SOLUTIONS 

The  numerical  solutions  to  the  five  experimental 
impact  situations  described  in  Table  I  are  presented  in  this 
section.  The  numerical  techniques  are  discussed  and  the 
material  model  employed  is  described  in  detail.  The  numeri¬ 
cal  results  are  then  presented  and  compared  to  determine  the 
effects  of  projectile  shape. 

4.1  NUMERICAL  TECHNIQUES 

The  HELP  code^  was  employed  to  solve  the  five  impact 
problems  which  make  up  the  projectile  shape  effects  study. 

HELP  is  a  multi-material,  two-dimensional,  advanced  Eulerian 
code  for  describing  material  response  to  dynamic  loading. 

Since  it  is  described  in  detail  in  Ref.  3,  only  a  brief  general 
outline  of  the  code  is  presented  here.  In  addition,  those 
features  of  HELP  which  are  specifically  pertinent  to  the 
current  study  such  as  interface  and  free  surface  treatment, 
closing  of  voids  and  the  inclusion  of  tracer  particles  which 
passively  trace  out  the  flow  field  history  are  discussed  in 
more  detail. 

Figure  6  shows  a  section  of  a  typical  Eulerian  grid 
used  in  the  HELP  calculations.  Since  the  y-axis  is  an  axis 
of  symmetry,  each  cell  actually  represents  the  cross  section 
of  an  annulus.  The  i-  and  j -lines  form  the  vertical  and 
horizontal  cell  boundaries  respectively.  The  governing  equa¬ 
tions,  i.e.,  conservation  of  mass,  energy  and  momentum  to¬ 
gether  with  constitutive  relations,  are  solved  throughout 

the  calculational  mesh  at  time  t.  The  process  involves 

r  3) 

three  phases  and  is  well  documented  elsewhere.  '  A  new 
time  increment,  At,  is  then  calculated  and  the  same  steps  are 
repeated  at  time  t  ♦  At.  Thus  the  solution  progresses 
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Fig.  6--Section  of  a  typical  Eulerian  grid  used 
in  the  HELP  calculations. 
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incrementally  in  time  and  the  material  response  history  can 
therefore  be  established. 

Special  features  have  been  incorporated  into  HELP  to 
minimize  certain  inherent  disadvantages  of  the  general  Eulerian 
technique  while  retaining  the  distinct  Eulerian  advantage  of 
handling  large  scale  material  flow.  For  example,  the  grid  can 
be  seeded  with  passive  tracer  particles  which  move  with  a 
weighted  velocity  of  the  material  in  their  vicinity.  This 
allows  a  particular  mass  element,  along  with  its  associated 
dynamic  variables,  to  be  traced  as  it  moves  through  the  cal- 
culational  mesh.  The  point  g  of  Fig.  6  represents  such  a 
tracer  particle.  The  dynamic  variables  associated  with  point 
g  are  calculated  from  weighted  averages  of  those  variables 
in  the  surrounding  cells  which  are  overlapped  by  the  ficticious 
tracer  particle  cell  shown  dashed  in  Fig.  6.  If  desired,  the 
tracer  particles  can  be  used  as  the  nodes  of  a  Lagrangian  type 
net  that  deforms  with  the  material. 

Interface  tracer  particles  such  as  point  "a"  of  Fig.  6 
are  employed  in  the  definition  and  propagation  of  material 
interfaces  through  the  Eulerian  grid.  This  represents  a 
distinct  advantage  over  the  general  Eulerian  technique  in 
that  material  diffusion  at  interfaces  is  not  allowed.  The 
special  interface  treatment  is  accomplished  by  calculating 
the  mass  flux  across  mixed  cell  boundaries  for  each  material 
and  imposing  the  restriction  of  like  material  flow  across 
partial  cell  boundaires.  The  interface  is  then  repositioned 
by  moving  the  interface  tracer  particles  according  to  a 
velocity  weighting  scheme  which  looks  across  the  interface. 

A  similar  scheme  is  used  for  free  surface  tracer 
particles  such  as  points  e  and  f  of  Fig.  6.  For  this 
case,  however,  a  portion  of  the  cell  is  void  so  that  only 
the  material  side  of  the  interface  influences  the  free  sur¬ 
face  velocity  weighting  scheme.  This  scheme  prevents 
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free  surface  diffusion  and  therefore  eliminates  another  in¬ 
herent  disadvantage  of  the  general  Eulerian  technique. 

For  the  projectile  shape  effects  study  a  special  auto¬ 
matic  void  closing  routine  was  incorporated  into  HELP.  This 
was  necessary  because  the  nature  of  the  pointed  projectile 
penetration  process  causes  initially  free  surface  regions  of 
both  the  target  and  projectiles  to  eventually  become  contact 
surfaces.  Hence,  two  free  surface  tracer  particles  such  as 
point  c  and  d  of  Fig.  6  may  later  become  a  single  inter¬ 
face  tracer  particle  such  as  point  b.  This  scheme  is  rather 
simple  in  that  it  mainly  involves  a  change  in  the  velocity 
weighting  scheme. 

The  foregoing  discussion  emphasizes  the  advantages  of 
HELP  over  the  general  Eulerian  method.  In  addition,  it 
describes  very  briefly  and  in  a  qualitative  manner  the  general 
numerical  scheme  as  well  as  the  special  features  employed  to 
handle  free  surfaces  and  material  interfaces.  For  a  very 
comprehensive  documentation  of  the  HELP  code  see  Ref.  3. 


4 . 2  MATERIAL  MODEL 


The  material  model  to  be  used  in  the  calculations 
consists  of  the  equation  of  state,  an  expression  for  the 
deviatoric  stress  components,  a  yield  criterion  and  a 
failure  criterion. 


The  equation  of  state  to  be  used  is  that  due  to  J.  H. 

Tillotson, modified  to  give  a  smooth  transition  between 

condensed  and  expanded  states.  For  the  condensed  states, 

i.e.,  when  p/p  >  1,  or  for  any  cold  states,  I  <  I  ,  the 

o  * 

equation  has  the  form 


Pc  * 


a  ♦ 


I 


+ 1 


Ip  +  Ap  +  By: 


(12) 
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For  expanded  hot  states,  i.e.,  when  p/p  <  1  and  I  >  I' 
the  equation  of  state  has  the  form 


PE  -  alp  + 


blp 


L1  n 


♦  l 


♦  Aye 


B(v/v  -1) 


-o(v/v  -1) 


(13) 


A  smooth  transition  between  the  condensed  and  expanded  states 
is  insured  by  a  transition  equation  for  the  intermediate  region 
defined  by  I_  <  I  <  I*  and  p/p  <  1.  This  blended  portion 
of  the  equation  of  state  has  the  form 


p .  (^Jpe  *  (£0 Pc 


(14) 


In  Eqs.  (12)  through  (14)  p,  I  and  p  are  pressure,  specific 

internal  energy  and  mass  density  respectively,  n  ■  p/p  *  y  ♦  1 

o 

-  v  /v,  and  p  ,  a,  b,  I  ,  I  ,  I\  A,  B,  o  and  6  are  con- 

0  0  0  5  5 

stants  for  the  particular  material.  The  values  of  these  con¬ 
stants  for  steel  and  aluminum  are  given  in  Table  II. 

The  deviatoric  stress  increments,  dS?.,  are  determined 

ij 

by  using  the  elastic  relation 

dS^  -  2Gdef.  (15) 

where  G  is  the  modulus  of  rigidity  and  de£j  are  the  incre¬ 
ments  of  deviatoric  strain.  When  such  an  increment  of  stress 
causes  the  yield  criterion  to  be  violated  each  stress  com¬ 
ponent  is  proportionately  reduced  to  bring  the  stress  state 
normally  back  to  the  yield  surface.  A  variable  yield  strength 

Y  .(Kf  ♦  Ktv  ♦  Kjti’)(l  -  I/I„)  (16) 
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can  be  defined  to  account  for  the  increase  in  strength  at 
high  pressures  and  the  decrease  in  strength  at  elevated 
temperatures . 

The  failure  criterion  to  be  employed  is  a  simple  one 
based  upon  relative  volume.  When  the  relative  volume  in  a 
cell  reaches  a  certain  value  which  is  greater  than  a  speci¬ 
fied  maximum  relative  volume,  that  cell  is  said  to  fail  and 
the  stresses  are  zeroed  out.  Table  III  gives  the  shear 
moduli,  yield  stresses  and  relative  volumes  for  failure  which 
were  used  for  aluminum  and  steel. 

4.3  NUMERICAL  RESULTS 


The  five  experimental  impact  situations  described  in 
Table  I  form  the  basis  for  the  numerical  projectile  shape 
effects  study.  The  corresponding  calculational  matrix  is 
given  in  Table  IV.  Since  the  projectile  nose  shape  was  of 
primary  interest  in  the  study,  materials,  target  thicknesses 
and  maximum  projectile  radii  were  not  varied.  The  initial 
projectile  kinetic  energies  are  not  identical,  however,  since 
it  was  of  interest  to  duplicate  the  experimental  conditions. 
The  last  column  of  Table  IV  indicates  the  number  of  usee  to 
which  each  solution  was  carried  out  in  time. 

Figure  7  is  a  composite  of  the  initial  projectile* 
target  configurations  employed  in  the  current  investigation. 
The  left-hand  side  of  each  projectile  represents  an  axis  of 
symmetry.  The  relative  dimensions  are  correctly  scaled. 

It  will  be  seen  later  that  this  method  of  presenting  the 
projectile-target  configurations  provides  a  convenient  means 
for  comparing  the  relative  penetration  depths,  regions  of 
plate  failure  and  projectile  deformations  as  the  solutions 
progress  in  time. 
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Since  many  of  the  general  qualitative  features  of  the 
five  numerical  solutions  are  similar,  only  the  solution  to 
the  Russian  projectile  impact  is  presented  in  detail.  The 
results  of  the  other  four  impact  problems  are  presented  in 
the  form  of  comparisons  of  specific  quantitative  data. 

Figure  8  shows  an  SC-4020  plot  of  the  initial  configu¬ 
ration  of  the  Russian  projectile  impact  problem  and  the 
Eulerian  grid  employed  in  the  solution.  For  clarity  in  sub¬ 
sequent  plots  the  grid  is  not  included.  These  plots  represent 
a  connection  of  the  free  surface  and  interface  tracer  particles 
and  hence  supply  the  projectile-target  configurational  history 
throughout  the  perforation  process. 

Figure  9  shows  the  configuration  at  2  usee.  The  shaded 
area  indicates  a  region  of  material  failure.  It  should  be 
noted  that  plate  failure  begins  along  the  axis  -*i  symmetry 
at  the  rear  free  surface  of  the  plate.  Figures  9  through  12 
indicate  that  the  failed  region  spreads  both  axially  and 
radially  indicating  petalling  failure.  A  complete  set  of 
the  configurational  plots  out  to  36  usee  is  provided  in  the 
Appendix. 

Because  it  would  be  impractical  to  include  complete 
sets  of  configurational  plots  for  all  five  solutions,  com¬ 
parisons  are  provided  in  Figs.  13,  14,  and  15  at  times  of 
10,  15  and  20  usee  respectively.  The  dashed  lines  in  these 
figures  indicate  the  initial  configuration  and  the  shaded 
areas  denote  regions  of  material  failure. 

In  Fig.  16  the  final  calculated  projectile  configura¬ 
tions  are  compared  with  their  undeformed  shapes.  The  times 
corresponding  to  each  of  these  configurations  are  given  in 
Table  IV. 
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Fig.  8- -Initial  configuration  and  computational  grid  for  the 
Russian  projectile  impact  calculation. 
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■MUM  (CM.  I 


Fig.  9--Russian  projectile-target  configuration  at  2  Msec. 
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*MIM  ICN.I 

Fig.  10--Russian  projectile-target  configuration  at  3  usee. 
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Fig.  ll>-Russian  proj ectile-target  configuration  at  5  usee. 
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V.  DISCUSSION  OF  RESULTS 

It  is  of  interest  to  compare  the  numerical  results 
with  both  perforation  foraula  predictions  and  experimental 
observations.  Figure  17  is  a  plot  of  normalized  projectile 
kinetic  energy  versus  time.  The  experimentally  observed 
residual  kinetic  eneigy  of  the  projectile.  Eg,  has  been  sub* 
tracted  from  the  numerically  predicted  projectile  kinetic 
energy,  (KE)^.  This  difference  has  been  divided  by  the  ex* 
pci. mentally  observed  loss  in  projectile  kinetic  energy, 

E0*Eg,  to  obtain  the  ordinate  for  the  plot  of  Fig.  17. 

Plotted  versus  time,  this  quantity  should,  for  all  projectile 
shapes,  have  a  value  of  unity  at  t  •  0  and  approach  zero 
asymptotically  et  some  later  time,  provided  the  HELP  predic* 
tions  are  in  agreement  with  experimental  observations.  This 
indeed  see as  to  be  the  case,  as  is  seen  in  Fig.  17.  The 
cylindrical  and  conical  nose  projectiles  attained  the  asympto* 
tic  values  rather  abruptly,  due  to  thair  bluntness.  The  more 
pointed  the  projectile,  the  longer  it  takes  to  attain  its 
correct  value.  This  is  because  the  projectile  must  travel  a 
greater  dista^c*  ]  if«re  the  maximum  projectile  radius  is  in 
contact  with  plate.  The  Russian  projectile  calculation 
uaj  carried  out  v  36  usee  before  it  approached  its  asymptotic 
value.  It  became  apparent  that  the  calculations  involving  the 
extremely  pointed  ogive  and  optimal  projectiles,  would  there* 
fore  have  to  br  carried  out  in  time  considerably  further 
than  36  usee.  Because  of  the  additional  computing  time 
necessary  to  carry  out  such  long  calculations,  it  was  decided 
to  stop  the  ogive  and  optimal  projectile  calculations  at 
29  wscc.  Based  on  the  fact  that  the  numerically  calculated 
cylindrical*,  conical-,  and  Russian-shaped  projectile  kinetic 
energies  all  attained  their  respective  experimentally  obser¬ 
ved  values,  it  is  believed  that  this  would  also  be  the  case  if 
the  ogive  and  optiml  impact  calculations  were  continued. 
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The  horizontal  lines  on  Fig.  17  labeled  T  and  N 
indicate  the  rhomson  and  Nishiwaki  theory  predictions.  If 
the  perforation  formula  predictions  were  in  complete  agree¬ 
ment  with  the  experimental  observations,  the  lines  would  all 
lie  on  the  time  axis,  i.c.,  they  *ould  have  zero  ordinate 
values.  Thus,  Fig.  17  indicates  that  the  numerical  predic¬ 
tions  o  residual  kinetic  energy  are  in  excellent  agreement 
with  experiment  but  that  the  Thomson  and  Nishiwaki  theories 
show  considerable  discrepancy.  The  values  of  T  and  N 
shown  in  Fig.  17  were  calculated  by  Fields^  using  the 
static  yield  strength  of  the  plate  material.  The  bracketed 
values  were  calculated  by  Good  using  the  Thomson  theory  and 
for  the  case  of  the  conical  projectile  there  is  considerable 
disagreement  with  Fields'  value. 

In  Fig.  18  the  percent  loss  in  projectile  kinetic 
energy  is  plotted  for  the  various  projectile  shapes,  beginning 
with  the  most  efficient  penetrator.  Th*  solid  circles  re¬ 
present  the  experimentally  observed  values.  The  open  circles 
represent  the  HELP  predictions.  As  expected  from  the  data 
presented  in  Fig.  17,  the  percent  loss  in  projectile  kinetic 
energy,  as  predicted  by  the  HELP  calculations,  for  the 
cylindrical  and  relatively  blunt  conical  head  projectiles 
are  in  good  agreement  with  experiment  while  the  more  pointed 
projectiles  are  not.  This,  of  course,  is  due  to  the  fact 
that  the  more  pointed  projectile  impact  solutions  were  not 
carried  out  sufficiently  far  in  time.  It  should  be  mentioned 
that  conversion  of  projectile  kinetic  energy  into  internal 
energy  was  still  taking  place  in  the  pointed  projectiles 
when  the  calculations  were  stopped.  Since  the  loss  in  pro¬ 
jectile  kinetic  energy  can  go  into  projectile  internal 
energy  as  well  as  target  kinetic  or  internal  energy,  energy 
partitioning  is  of  interest.  Figure  18  shows  the  numerically 
predicted  energy  partitioning  by  means  of  shaded  regions. 
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Fig.  18--Percent  loss  in  projectile  kinetic  energy  for  the 
various  projectile  shapes.  The  arrows  attached  to  the  open 
circles  indicate  that  the  calculated  projectile  kinetic 
energies  were  still  decreasing  when  the  calculations  were 
terainated. 
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The  broken  lines  connecting  the  open  triangles  and 
squares  in  Fig.  18  show  the  percent  reduction  of  projectile 
kinetic  energy  predicted  by  the  Thomson  theory  using  the 
dynamic  and  static  target  yield  strengths  respectively. 

The  predictions  of  the  Nishiwaki  theory  are  shown  for  the 
cylindrical  and  conical  projectile  impacts  by  the  dashed 
line  connecting  the  solid  triangles. 

The  solid  square  shows  the  value  predicted  by  the 
Thomson  theory  with  static  yield  strength  as  reported  by 
Fields^  for  the  conical  head  projectile  impact.  The  value 
of  percent  reduction  in  kinetic  energy  obtained  by  Good  who 
applied  the  same  formula  was  considerably  greater  as  evidenced 
in  Fig.  18.  Since  Fields  did  not  calculate  energy  loss  for 
the  Russian  projectile,  his  value  for  the  conical  projectile 
lead  him  to  state  that  “Preliminary  calculations  using  the 
Thomson  theory  predicted  very  small  differences  in  perfora¬ 
tion  velocity  loss  between  the  different  shapes."  ^  The 
dashed  line  which  connects  the  open  squares  in  Fig.  18  shows 
that  the  percent  loss  in  projectile  kinetic  energy  predicted 
by  the  Thomson  theory  for  both  the  Russian  and  conical  pro¬ 
jectile  shapes  indicate  disagreement  with  Fields'  statement. 
It  is  concluded  that  the  Thomson  theory  is  sensitive  to  pro¬ 
jectile  shape  for  the  four  pointed  impact  cases  reported 
here,  but  that  the  predicted  values  of  projectile  energy  loss 
are  low. 

The  residual  velocity  for  the  cylindrical  impact 
case  calculated  by  Heyda's  formula  is  4.43  *  10*  cm/sec  as 
compared  to  4.61  «  10*  cm/sec  predicted  by  Nishiwaki 's  formula 
and  4.S  *  10*  cm/sec  predicted  by  Recht's  formula  using  a 
ballistic  limit  of  1.02  «  10*  cm/sec.  The  experimentally 
observed  value  of  residual  velocity  of  the  cylinder  was 
4 . 54  *  10*  cm/sec.  Although  Recht's  formula  seems  to  be 
closest  to  the  observed  value  for  this  case,  the  total  spread 
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in  predicted  velocities  is  less  than  4  percent  of  the 
observed  value. 

Figure  18  also  provides  a  aeasure  of  relative  pene¬ 
tration  efficiency  for  the  various  projectile  shapes.  Of 
the  five  shapes  investigated,  the  optiaal-head  projectile 
is  the  most  efficient,  followed  by  the  ogival*,  Russian*, 
cylindrical-  and  conical-head  projectiles. 
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VI.  CONCLUSIONS 

Within  the  scope  of  the  current  investigation, 
available  pointed  projectile  perforation  formulas  do  not 
adequately  predict  residual  projectile  kinetic  energy. 

HELP  calculations,  carried  out  sufficiently  far  in  time, 
seem  to  predict  values  of  projectile  kinetic  energy  which 
asymptotically  approach  the  observed  residual  energies. 
This  agreement  with  experimental  observations  is  quite 
encouraging,  since,  in  addition  to  adding  physical  insight 
into  penetration  phenomena,  the  results  demonstrate  that 
numerical  solutions  can  provide  useful  information  to  both 
weapons  and  armor  designers,  as  well  as  to  vulnerability 
analysts.  Furthermore,  the  HELP  code  should  be  a  useful 
tool  in  parametric  studies  designed  especially  for  de¬ 
veloping  new,  or  altering  presently  available,  simple 
perforation  formulas. 
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RUSSIAN  PROJECT I LE -TARGET  CONFIGURATIONS 
AT  VARIOUS  TIMES 
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Fig.  A.2-*Russian  projectile  tirgct  configuration 
at  t  •  2  nsec. 
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Fig.  A.3--Russian  projectile- target  configuration 

at  t  •  3  uaec. 
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Fig.  A.4--Russian  projactila-targct  configuration 

at  t  -  5  ysoc. 
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Pig.  A.5--Ru»*itn  projactile-targat  configuration 

at  t  ■  10  ysac. 
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Fig.  A 


.b-'Russitn  proji ctilc-ttrgit  configuration 
at  t  ■  15  usee. 
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Fig.  A.7--Russian  projoctile-target  configuration 

at  t  ■  20  u*ec. 
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Fig.  A.8--Russian  projactila-targat  configuration 

at  t  "  30  ysac. 


Fig.  A.9--Russian  projectile-target  configuration 

it  t  ■  36  usee. 


